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ABSTRACT 

Co-addition of deep (rms ~ 30/iJy) 20cm data obtained with the Australia Tele- 
scope Compact Array at the location of Spitzer Wide field survey (SWIRE) sources 
has yielded statistics of radio source counterparts to faint 24/xm sources in stacked 
images with rms < 1/iJy. We confirm that the infrared-radio correlation extends 
to /24pm = 100/iJy but with a significantly lower coefficient, /20cm — 0.039/24/im 
(924 = log(/24^m//20cm) = 1.39 ±0.02) than hitherto reported. We postulate that this 
may be due to a change in the mean 924 value ratio for objects with /24/xm < lmJy. 
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1 INTRODUCTION 

The tight correlation between the global far-infrared and 
radio emission from galaxies has been established for over 
three decades. Early ground-based studies (van der Kruit 
1973; Condon et al. 1982) and observations with the infrared 
Astronomical Satellite (IRAS, Dickey & Salpeter 1984; de 
Jong et al. 1985) conclusively established a correlation over 
a broad range of Hubble types and luminosities at low red- 
shifts (z < 0.1). The origin of this correlation is thought 
to lie in the link between massive stars, which generate far- 
infrared emission by reheating dust, and supernovae, which 
accelerate cosmic rays that then generate radio synchrotron 
radiation (Harwit & Pacini 1975; Condon, 1992; Xu, Lisen- 
field & Volk 1994). However, alternative mechanisms have 
also been proposed, including secondary electron production 
(radio synchrotron) from molecular clouds (far infrared). 

More recently, by extrapolating observations made with the 
Infrared Space Observatory (ISO) in the mid-infrared (Co- 
hen et al. 2000, Garrett 2002) and in the sub- millimetre with 
the SCUBA instrument, (Carilli & Yun 1999, Ivison et al. 
2002), it has been argued that the far-infrared/radio corre- 
lation extends to z > 1. If confirmed, this suggests that the 
galaxies at high redshift may share many of the same prop- 
erties as their low redshift counterparts; whatever the origin 
of this relation. 

With the greatly enhanced mid-infrared sensitivity provided 
by the Spitzer Satellite, it is now possible to investigate di- 
rectly the far- infrared/radio correlation at redshifts of order 
unity. An initial attempt has already been made by Ap- 
pleton et al. (2004), using the Spitzer First Look Survey 
(FLS) and the Very Large Array (VLA). The flux limits 
for the FLS and VLA samples were /24 M m = 0.3mJy and 



/20cm = 115mJy respectively. Although limited by the avail- 
ability of spectroscopic redshifts, Appleton et al. were able 
to demonstrate a good correlation between the 24/im and 
20cm flux for objects with a median 2 = 0.3 extending out to 
z ~ 1. Their analysis was limited by the relative brightness 
of the radio flux limit, which was a factor of three brighter 
than that corresponding to the 24/^m limit, based on the 
observed far-infrared/radio flux ratio. 

In this analysis, we extend the flux limits to which the far- 
infrared/radio correlation may be studied by a further factor 
of three at 24pm and by a factor of 20 at 20 cm, using the 
Spitzer Wide Field Survey (SWIRE) and the technique of 
stacking (see e.g. Wals et al. 2005) applied to deep Australia 
Telescope Compact Array (ATCA) observations. 

2 DATA 

2.1 ATCA observations 

Two regions around the Chandra Deep Field South (CDFS, 
RA(2000)=3 h 31 m Dec(2000)=-28°06 m ) and the European 
Large Area ISO Survey SI (ELAIS, RA(2000)=0h34, 
DEC(2000)=-43°45 m , have been observed with a series of 
mosaiced pointings at 20cm with the ATCA as part of the 
Australia Telescope Large Area Survey (ATLAS, Norris et 
al. 2006; Middelberg et al. 2007, in prep). The goal of the 
ATLAS program is to provide large area (6 deg 2 ), deep sur- 
veys of radio sources in fields which are also the subject of 
intensive study at other wavebands. In particular, the areas 
covered are chosen to be matched to those covered by the 
SWIRE survey (Lonsdale et al. 2003) as well as enabling op- 
tical spectroscopy with instruments such as the Anglo Aus- 
tralian Telescope 2-degree field. A total of 360 hours inte- 
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Figure 1. ATCA 20cm images of the CDFS field (left) and the ELAIS field (right). The areas used in these analyses is denoted by the 
circles. 



gration has been obtained on the CDFS, and of 234 hours in 
the ELAIS. In each field, 28 and 20 pointings (each with a 
mean integration of 13 hours) have been mosaiced together 
to provide an area of approximately 1.6 deg x 2.2 deg in the 
CDFS, and of 1.6 deg x 1.6 deg in the ELAIS. 
The data were calibrated and imaged using Miriad (Sault 
et al. 1995) following the procedures in the Miriad cook- 
book for high dynamic range observations. In particular, 
mult i- frequency clean had to be used to account for the 
high fractional bandwidth of 15 per cent. The mean rms 
noise in the images is 30/iJy, and is limited not by integra- 
tion time but by sidelobes arising from strong sources in 
(CDFS) or just outside (ELAIS) the imaged areas (shown 
in Fig 1). The CDFS image also contains data from an 
earlier observation by Koekemoer et al. (in preparation), 
who have covered a region of 1 deg diameter centred on 
RA(2000)=3 h 33, Dec(2000)=-27° with a mosaic of seven 
pointings, and reach an rms noise level of ~ 15/iJy. In or- 
der to avoid regions around the edge of the mosaic with 
increased noise levels, we restricted the analysis to a 48- 
arcmin radius centred on the ATLAS field centres, denoted 
with circles in Fig 1. As a further check of the calibration, 
we confirmed that there was no significant offset between 
the published NVSS fluxes and the flux estimates obtained 
from the ATCA data for bright sources (/20cm > 2 mjy) in 
the CDFS field (see also Norris et al. 2006) 



2.2 Spitzer observations 

The Spitzer Space Telescope was launched in August 2003. It 
is equipped with several instruments, and here we use data 
taken with the Multiband Imaging Photometer for Spitzer 
(MIPS) in its 24/im band. Before launch, proposals were 
invited for Legacy Science Programs, which were large co- 
herent programs whose data would be of lasting importance 
to the broad astronomical community. One of the six Legacy 
projects chosen was the SWIRE (Spitzer Wide Extragalac- 
tic) Survey (Lonsdale et al. 2003), which has observed a 
region of 6 square degrees surrounding the CDFS, and of 5 
square degrees coincident with the ELAIS field. The analy- 
sis of those data are described by Surace et al. (2006) . Most 
of the SWIRE data are now in the public domain (SWIRE 



data Release 3 1 ). Here we use only the 24/nm MIPS data. 
The fluxes used here are, in the case of unresolved or noisy 
sources, aperture-corrected fluxes, as described by Surace et 
al. (2006). In the case of extended sources, Kron fluxes are 
used (Kron 1980). 

3 STACKING ANALYSIS 
3.1 Method 

We present here in detail the analysis of the radio and in- 
frared data obtained in the CDFS and ELAIS field. Both 
fields were analysed separately. A more detailed analysis was 
conducted on the CDFS field because of the more extensive 
data set (deeper SWIRE data, optical CCD data), but the 
ELAIS field provided a useful independent check of the re- 
sults. 

We first removed galactic stars from the SWIRE catalogue 
in both fields. For the CDFS we used the optical-infrared 
colour (r — i/r — 3.6/im) selection method employed by 
Rowan- Robinson et al. (2004). Optical magnitudes for all 
SWIRE sources were obtained from CTIO CCD observa- 
tions (Siana et al. in preparation). At bright optical magni- 
tudes (r < 17), image saturation renders the optical magni- 
tudes unreliable and so the UKST R-band sky survey images 
of bright SWIRE (/24 M m > 1000/iJy) sources were inspected 
visually using the SuperCOSMOS sky server to remove any 
further stellar contamination. For the ELAIS field we did 
not have access to reliable broadband optical magnitudes 
and so relied on the stellar classification provided by the 
SuperCOSMOS catalogue (confirmed by visual inspection 
for the brighter images). Reliable star/galaxy classification 
is only possible at r < 20mag from these measurements, 
and so it is possible that some residual stellar classification 
remains amongst the fainter optical counterparts. 
Using AIPSH — h we then extracted 43 x 43 arcsec 2 sub- 
images from the 20 cm images, centred on each remaining 
SWIRE extragalactic source above the 5u 24/im flux limit in 
the CDFS (/ 2 4 M m > lOO^iJy) and ELAIS (/ 2 4 M m > 400/iJy) 
fields contained within the 48-arcmin radius region. 
The 20cm images were then sorted into logarithmic (0.1 dex) 

1 http://swire.ipac.caltech.edu/swire/astronomers/data_access.html 
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Table 1. Flux measurements of stacked images in the CDFS field 
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24^im flux bins. We computed a stacked radio image for each 
bin by forming the median of corresponding pixels in all ra- 
dio images in that bin. For each bin, two stacks were cre- 
ated: one which included all sources (the all-source stack) 
and one which excluded those images whose central pixels 
had /20cm > 100/xJy (the quiet-source stack). The radio flux 
for each bin was then estimated from the measured peak flux 
in the stacked image. We also measured fluxes by fitting a 
point source to the image. This yielded identical results, in- 
dicating that the stacked images were not resolved. We also 
confirmed that stacking similar number of random positions 
within CDFS field produced no significant detections. Im- 
ages of the quiet-source stacks are shown in Fig 2. 
To estimate the accuracy of the flux determinations, we 
created simulated radio images in AIPS+- 1- with constant 
noise levels (rms=30^iJy) across the CDFS field. To this we 
added sources at the positions of all SWIRE sources with 
/24 M m > lOO^iJy with a scaled, / 2 o C m = 0.04/ 2 4 Mm , radio 
flux. We repeated the analysis above, computing both all- 
source and quiet-source stacks. The simulated median im- 
ages are shown in Fig 3. As a further check on our analysis 
processes, we repeated the simulation by adding the artifi- 
cial sources to simulated u-v data, which we then processed 
in an identical way to the real u-v data. This yielded an 
essentially identical result to the simulations that had been 
performed in the image plane. 



3.2 CDFS field 

The resultant fluxes for the stacked images in the CDFS 
are listed in Table 1, together with the median infrared 
source flux /24 Mm in each bin. In many of the fainter /24 Mm 
bins, we have been able to stack over 1000 images. How- 
ever, we note that the rms in the resultant stacked image 
(estimated from the image by excluding the central quarter) 
has not decreased in line with Poisson statistics, probably 
because of confusion from sub-^ijy radio sources. Assuming 
the infrared/radio correlation identified below, sources with 
/24pm = 10/iJy will have /20cm ~ 0.5piJy, and extrapolating 
the /24/jm number counts (see Table 1) predicts a source den- 



sity of over 10 sources deg~ at /24/jm = 10/xJy, approaching 
10 sources per stacked image. A confusion noise threshold 
of ~ 0.5/iJy when added in quadrature to an assumed Pois- 
sonian reduction in the noise afforded by the stacking pro- 
cedure, i.e. 30/xJy / 'y/n, is consistent with the observed rms 
noise measurements. 

The images from the quiet-source stacks in the CDFS field 
are shown in Fig 2, and the median values of f2i^.m and /20cm 
are plotted in Fig. 4. Corresponding images of simulated 
data are shown in Fig 3, and their /20cm vs /24pm plots 
are shown in Fig. 5. The qualitative agreement between the 
corresponding real and simulated images is good but not 
excellent. The simulations do not reproduce the fine scale 
noise in the real data, which we believe is partially due to 
the residual effects of the bright source to the southwest. In 
support of this conclusion, we note that by moving the field 
centre closer to that source the level of the fine scale noise 
increases slightly in the stacked images. 
Nevertheless, we note that the flux estimates from the sim- 
ulated data (both the all-source and quiet-source stacks) are 
fitted by a linear relation 

/20cm = 0.044/ 24M m + 4.74^Jy 

which is essentially identical to the input simulation. We 
conclude that the stacking introduces no systematic bias 
into our flux measurements. 

We note that the /20cm v /24pm correlation flattens off at 
the brightest flux bins for the quiet-source stack. This re- 
veals the systematic under-estimation of the radio flux in 
the quiet-source stack caused by the exclusion of the increas- 
ingly large number of radio sources (> 20 per cent) at the 
/20cm > 100/iJy level that are bona-fide counterparts to the 
SWIRE sources at high /24pm fluxes. Note that this bias 
is also seen for the simulated quiet-source stacks, although 
only significant for simulated stacks with /24 Mm > 2000 /iJy. 
We obtain a fit to both the all-source stack (all bins) of the 
form: 

/20cm = 0.041(±0.002)/ 24M m + 1.35(±0.8)/xJy 
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Figure 2. Images from the quiet-source stacks. Upper left bin 
corresponds to faintest bin; 2 < log/24 (Jm < 2.1. Bins increase 
log/24/jm =0.1 from left to right and top to bottom. 



Figure 3. Images from the simulated quiet-source stacks. Upper 
left bin corresponds to faintest bin; 2 < log/24 Mm < 2.1. Bins 
increase log/24/jm = 0.1 from left to right and top to bottom. 



or an infrared-to-radio flux ratio, ^24 = log(/24,im//20cm) = 
1.39 ± 0.02. The data points derived for quiet-source stack 
have a slightly lower normalisation throughout the infrared 
flux range. This is because objects are detected above our 
radio limit in all /24 M m bins used in this analysis. The frac- 
tion with radio detections rises from 2 per cent for the lowest 
/24 M m values to 40 per cent at the highest (see Table 1). 



3.3 ELAIS field 

Following the same procedures we obtained a value for the 
infrared- radio correlation using the ATCA and SWIRE ob- 
servations of the ELAIS region. The ELAIS region does not 
extend to as deep a 24/itm flux limit as the CDFS SWIRE 
field (see Table 2). A fit to the all-source stack gives the 
relation: 



/20cm = 0.039(±0.004)/ 24Mm +7.1(±3.3)MJy 

for the all-source stack. The slope is consistent with the 
result obtained from the CDFS data, giving a §24 = 1.41 ± 
0.04. The intercept is different but this is largely dominated 
by fitting to the data points at the lowest flux levels where 
the errors are high and the fit poorly constrained. Forcing 
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Figure 4. The correlation between median radio flux, /20cm, and 
infrared flux, /24pm, for the stacked CDFS sources. Flux mea- 
surements are included for stacks with (filled circles) and without 
(open circles) the inclusion of detected (> lOO^tJy) radio flux at 
the SWIRE sources positions. The numbers denote the numbers 
of SWIRE source positions included in each stack). The dashed 
line is the best- fit line to the data obtained from the 'all-source' 
stack. 
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Table 2. Flux measurements of stacked images in the ELAIS field 
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Figure 5. The correlation between median radio flux, /20cm? and 
infrared flux, /24^m, for the simulated data with the scaling 0.04. 
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Figure 7. A composite plot of the CDFS field shown /24,zm and 
/20cm fluxes for individual sources (small triangles) and the all- 
source (filled circles) and radio-quiet (open circles) stacked im- 
ages. 



the correlation to go through the origin of the /24 Mm / /20cm 
relationship gives an consistent fit in both cases with 924 = 
1.40 ±0.02. The correlation is shown in Fig. 6 and the fluxes 
for the stacked images are given in Table 2. 



100 . 1 000 

F 24,m <™ 

Figure 6. The correlation between median radio flux, /20cm, and 
infrared flux, /24/am, for the stacked ELAIS sources. Flux mea- 
surements are included for stacks with (filled circles) and without 
(open circles) the inclusion of detected (> 100/^Jy) radio flux at 
the SWIRE sources positions. The numbers denote the numbers 
of SWIRE source positions included in each stack). The dashed 
line is the best-fit line to the data obtained from the 'all-source' 
stack. 



4 DISCUSSION 

The most significant result to emerge from this analysis is 
that infrared-selected sources broadly follow the radio/far- 
infrared (FIR) correlation down to microjy flux densities, 
implying that these sources are primarily driven by star- 
formation rather than by AGN activity. This result rules out 
models such as that of Wall et al. (1986), who propose that 
the low-flux radio sources are low-luminosity radio galax- 
ies. Moreover, it confirms predictions that the far-infrared 
emitting starbursts should dominate the radio population at 
/20cm < 50/iJy (Gruppioni et al. 2003) based on extrapola- 
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Figure 10. MIPS 24/im images for four CDFS source with 
0.8 < /24pm < 0.83mjy and no detected 20cm emission. Note the 
complex/extended morphology in the majority of cases, broadly 
consistent with the population at this flux level and 924 value. 



Figure 8. Upper: Radio flux, /20cm: and infrared flux, /24um, for 
the individual CDFS sources. Fits using Lower: As for upper plot, 
with simulated sources: half with /20cm = 0; half with /20cm = 
0.04/24 M m, based on CDFS /24^m distribution. See text for more 
details. 



tions of models which fit the mid-infrared/radio correlation 
at brighter fluxes. 

The results presented here show that, in broad terms, this 
same correlation holds for sources down to microjy lev- 
els, thus extending the range over which the correlation is 
known to hold by some two orders of magnitude in flux den- 
sity. However, the data do not immediately tell us whether 
these galaxies are high-luminosity high-redshift star-forming 
galaxies or low-luminosity low-redshift star-forming galax- 
ies. 

Most previous studies of the radio-FIR correlation have 
measured qpir = log(S l FiH,/S , 20cm), where Sfir is an inte- 
grated far infrared flux typically obtained from IRAS fluxes 
using a weighting algorithm e.g. SpiR^Wnr 2 ) = 1.26 x 
10 4 (2.58S , 6i)Mm + <Sioo/jm)Jy. Such studies, whether on radio- 
selected (Condon et al. 1991), infrared-selected (Condon & 
Broderick 1991), or distance-limited (Helou et al. 1985) sam- 
ples, all broadly give the same value of qFiR = 2.3 ± 0.2. 
Here we measure q 2 4 = log (/24 M m //20cm), rather than qFiR, 
so to compare our results with previous authors requires a 
conversion factor between §24 and qFiR, which will depend 
sensitively on the SED of the parent population and the 
MIPS filter response, and is beyond the scope of this paper. 
However, 924 has been measured by Appleton et al. (2005) 
who obtain a (non-k-corrected) value of 524 = 0.84, which 
differs significantly from our measured value of ^24 = 1.39. 
Although our value of q is derived from stacked data, whilst 
the Appleton value is derived from detected sources, this 
disrepancy is nevertheless a surprising result and merits fur- 
ther scrutiny. First, we note that the simulations discussed 
in Section 3 have eliminated the possibility that it might be 
caused by an error in our imaging, stacking, or source-fitting 
algorithms or software. 

Another possibility is that the discrepancy might be caused 
by a difference in k-correction between our analysis and that 
of Appleton et al. (2004). We have made no k-correction to 
this flux ratio because we have neither redshifts nor detailed 



spectral information on these objects in either the infrared 
or radio bands. Appleton et al. (2004) explicitly demonstrate 
that individual k-corrections have little systematic effect on 
the mean trend of g 2 4 as a function of redshift at z < 1. 
Based on spectral energy distribution modelling in the in- 
frared, and an assumed synchrotron spectral index of —0.7 
in the radio, Appleton et al. found that this increased the 
uncorrected value of §24 = 0.84 by only 0.1-0.2, depending 
on the SED template used. 

However, Fadda et al. (2006) demonstrate that the k- 
correction arising from the PAH emission features enter- 
ing the 24/mi band at higher redshifts (1.1 < z < 1.8 or 
2.2 < z < 2.5) could be much greater. Fadda et al. (2006) 
found that emission in the 24pm band could be enhanced by 
as much as a factor 2 (5q = 0.3) for an M82 starburst galaxy 
template. Thus it is possible that k-corection provides a par- 
tial explanation for the differences seen with Appleton et 
al. (2004), if the bulk of the population used in the stack- 
ing analysis are high redshift starbursts. For k-correction to 
fully explain the difference our sources would need to be 
dominated by high redshift (z > 1) starburst galaxies with 
PAH emission features approximately double that seen in 
some of the most intense star-forming galaxies locally. 
A more mundane source for the discrepancy may arise from 
the different sample selection effects. We note that <?24 de- 
rived by Appleton et al. (2004) is based on radio detections 
of infrared sources. Our values derived here are for all in- 
frared sources, irrespective of radio flux. Thus a lower coef- 
ficient could be explained by a population of sources with 
low radio-infrared flux ratios which have previously not been 
included in the study of radio-selected samples because they 
fall below the radio flux limit for the sample. 
Indeed, we note that Gruppioni et al. (2003) found that the 
mid-infrared (15/Ltm) /radio (20 cm) luminosity correlation 
for mid-infrared sources with radio detections was a factor 
of two higher than the correlation for all sources with upper 
limits for radio luminosities included in the analysis. How- 
ever, if this were to be the case we would still have to ex- 
plain the difference between these lower 924 (passband issues 
notwithstanding) and those FIR/radio correltions derived 
from IRAS sources with complete radio detections (Yun, 
Reddy & Gruppioni 2003). 

To investigate this further, we plot in Fig. 7 a composite 
of the radio/infrared fluxes derived for both the individ- 
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Figure 9. ii-band UK Schmidt telescope images of nine galaxies with anomalously low radio-infrared flux ratios. Each image is 30x30 
arcsec. For top-left to bottom-right, the galaxies have /24^ m fluxes of 13.5, 6.99, 6.06, 5.97, 5.75, 5.29, 3.98, 3.97 and 3.86mjy, and all are 
below the 3<r radio flux limit (/20cm < 100/iJy). As indicated, four galaxies have spectroscopic redshifts from the 2dF Galaxy Rcdshift 
Survey (Colless ct al. 2001). 



ual radio-detected sources in the CDFS (Norris et al. 2006) 
and the stacked images derived in this paper. The sample 
of detected sources includes both AGN and star forming 
galaxies, resulting in a significant tail in Fig. 7 extending 
to high radio-FIR ratios, and precluding a measurement of 
q24- However there is a significant absence of sources below 
/20cm = 0.1/24 M m, (the radio-FIR relation found by Apple- 
ton et al.), which appears to mark a lower bound to the 
distribution. In particular, there is a significant absence of 
sources along the line extrapolated from our stacked data. 
Indeed, at the higher (/24 Mm > lOmJy) fluxes nearly all non- 
stellar 24/im sources are detected at 20cm in the CDFS field 
(/20cm > 0.4mJy). 

As a further consistency check of our result we derived the 
infrared-radio correlation based on radio and infrared fluxes 
measured at the positions of individual sources used in our 
stacking sample. Most of the radio sources will lie below 
the formal 5a detection limit and consequently their fluxes 
will be subject to large errors. Nevertheless the ensemble 
properties of the distribution will be robust. In each infrared 
bin used above we computed the median value of the fluxes 



in each bin and made the fit to those points. For the CDFS 
and ELAIS the resulting fits were respectively: 

/20cm = 0.037(±0.001)/24 M m + 2.6(±0.8)^iJy Median 

/20cm = 0.041(±0.002)/ 24M m + 6.6(±3.5)^Jy Median 

The median estimator to the individual points in both fields 
gave a consistent slope to that obtained from the stacked 
images. 

We also modelled the effects of a population with a bi- 
modal radio- infrared flux ratio; half of the population having 
/20cm = and half having /20cm = 0.08/24 M m (based on the 
CDFS MIPS flux distributions). Although the stacked im- 
ages of this simulated population give similar results to that 
obtained with the real data, the distribution of the individ- 
ual simulated population flux values is visibly bimodal in 
a scatter plot of the individual source fluxes (Fig. 8 upper 
panel) and clearly different from a similar plot using the real 
data (Fig. 8 lower panel) 

Nevertheless, we are faced with the conundrum that, what- 
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ever the distribution, the mean observed 524 value is greater 
than that of any modelled SED. Some of the optically- 
brightest CDFS sources with the most extreme 924 values 
are relatively low redshift galaxies (Fig. 9). However, most 
of these radio-weak sources have no direct optical counter- 
part at photographic plate limits (R ~ 21 mag). The major- 
ity of the objects at high /24pm flux end of the correlation 
/24pm ~ 0.83mJy and show extended or complex morpholo- 
gies in their 24pm images (Fig. 10). This would be consistent 
with a population dominated by actively star-forming galax- 
ies at moderate redshifts. To date, optical spectroscopy has 
been obtained for w 130 SWIRE sources in the CDFS (Nor- 
ris et al. 2007). The median redshift is 2 ~ 0.4 with a tail to 
higher redshift 2 ~ 2, exclusively composed of broad emis- 
sion line active galactic nuclei at z > 1.2. However, this red- 
shift distribution is strongly influenced by optical selection 
effects - the vast majority of the objects with measured red- 
shifts having relatively bright optical magnitudes (R < 22.5 
mag). 



5 CONCLUSIONS 

We find that the population of objects at /ijy flux densities 
obey a radio- 24pm correlation which indicates that they are 
powered by star formation rather than by AGN, but their 
mean value of 524 is inconsistent with previous results. We 
therefore propose that at microjy sensitivity, radio surveys 
are sampling a different population of objects, or objects 
with an intrinsically different value of 524, from objects seen 
in the local Universe. Nevertheless, we acknowledge that the 
change obtained here for the mean value of 524 occurs when 
we move from the use of discrete radio (> 5a) detections 
to fainter fluxes. Despite extensive testing and simulations 
we cannot rule out some unsuspected selection effect which 
may also contribute to this result. A partial explanation of 
the discrepancy may arise from the k-correction, but only 
if the population is dominated by extreme starburst galax- 
ies at 1 < 2 < 1.8 or 2.2 < z < 2.5. Another explana- 
tion for at least part of the discrepancy may arise from 
the radio-selection biases inherent in some previous surveys. 
Data from high-sensitivity surveys such as ATLAS will help 
to resolve such issues. 
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